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DIELECTRIC  RELAXATION  AND  ROCK 
GEOPHYSICAL  CHARACTERISTICS 


SECTION  I 

PROJECT  OBJECTIVES 


U 


i.j 


Li 
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The  ultimate  objective  is  to  adopt  rock  dielectric  parameters  as  indicators 
of  their  geophysical  characteristics  and  genesis.  To  achieve  this  goal, 
dielectric  relaxation  mechanisms  in  representative  rocks  must  be  developed. 
Electrical  models  simulating  rock  dielectric  dispersions  will  be  utilized  to 
correlate  rock  properties  with  the  model  parameters.  Self-  and  mutual- 
impedance  data  of  rocks  at  a  broad  frequency  range  will  be  gathered. 


RESEARCH  OBJECTIVES  FOR  THE  FIRST  SIX  MONTHS 

The  objectives  during  this  period  were  to  modify  our  existing  impedance 
equipment  with  its  on-line  computer  capabilities  to  measure  the  complex 
impedance  of  rocks  in  the  frequency  range  0.  05  to  2,  000  Hz.  The  next  step 
was  to  conduct  impedance  measurements  on  three  selected  rock  samples 
provided  by  the  contracting  agency  and  to  display  the  data  on  an  Argand 
Diagram.  A  further  objective  was  to  develop  techniques  for  the  transforma¬ 
tion  of  impedance  parameters  to  dielectric  parameters  and  to  construct 
Cole-Cole  plots  for  a  representative  rock. 


SIGNIFICANT  RESULTS 

An  analog/ digital  computer  setup  was  assembled  to  evaluate  rock  impedance 
in  the  frequency  range  0.05  to  2  kHz.  Room  temperature  neasurements 
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were  made  on  several  samples  of  basalt,  granite  and  quartz.  The  data 
appear  to  verify  the  impedance  circular -arc  diagram  when  displayed  in  the 
4rgand  plane.  Figures  1  through  3  show  the  actual  computer  output  for  the 
impedance  components  of  basalt,  granite  and  quartzite  rock  samples. 

The  adherence  of  the  data  to  the  circular-arc  plot  is  clearly  evident.  Quali¬ 
tative  effects  were  found  for  both  temperature  and  rock  humidity.  A  method 
was  developed  to  correct  the  rock  impedance  experimental  data  for  electrode 
polarization  artifacts.  A  new  technique  was  developed  to  extract  the  rock 
dielectric  constants  from  their  corrected  and  normalized  impedance  parameters. 
Detailed  results  and  their  analyses  are  given  in  Section  III  of  this  report.  The 
technical  report  summary  is  given  in  Section  IV. 
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Figure  1.  Computer  Display  of  the  Impedance  Circular 
Arc  for  a  Dresser  Basalt  Sample.  Points 
(+)  represent  experimental  variation  of 
reactance  with  resistance  at  a  given  frequency. 
Dotted  curve  is  an  arc  of  the  computer  data- 
filled  circle  in  the  least  square  sense. 
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Figure  2.  Computer  Display  of  the  Impedance  Circular  Arcs 
for  Granite  Samples.  Top  photograph  is  for  a 
long  granite  cylinder  (4.41  cm).  Bottom  photo¬ 
graph  is  for  a  thin  slice  of  the  granite  cylinder 
(0.  15  cm).  Points  (+)  represent  variation  of 
reactance  with  resistance  at  a  given  frequency. 
Dotted  curves  are  arcs  of  the  computer  data- 
t’itted  circle  in  the  least  square  sense. 
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Figure  3.  Computer  Display  of  the  Impedance  Circular  Arc 
for  a  Large  Diameter  Quartzite  Disc  Sample. 
Points  (+)  represent  experimental  variation  of 
reactance  with  resistance  at  a  given  frequency. 
Dotted  curve  is  an  arc  of  the  computer  data- 
fitted  circle  in  the  least  square  sense. 
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SECTION  II 

DESCRIPTION  OF  EXPERIMENTAL  EQUIPMENT  AND  TECHNIQUE 

A.  EXPERIMENTAL  SETUP 

The  basic  approach  to  the  experimental  determination  of  the  rock  impedance 
is  the  use  of  the  rock  as  a  feedback  impedance  in  an  operational  amplifier. 
Thus,  for  a  known  input  signal  and  resistance,  the  rock  impedance  can  be 
determined  from  the  amplifier  output  signal.  This  technique  lends  itself  to 
convenient  on-line  signal  analysis  through  the  use  of  analog- to -digital  con¬ 
version  equipment. 

Figure  4  illustrates  the  experimental  setup  in  simplified  form.  Basically, 
a  sinusoid  at  known  frequency  drives  the  amplifier  containing  the  rock,  while 
an  in-phase  square  wave  generated  simultaneously  with  the  sinusoid  is  used 
to  clearly  define  the  period  of  the  wave  form. 

The  output  from  the  rock  amplifier  is  sampled  over  some  interval  consistent 
with  the  signal  frequency  and  converted  to  digital  form.  At  this  state  a 
DDP-24  computer  performs  a  discrete  Fourier  transform  and  computes  the 
equivalent  series  resistance  and  reactance  for  the  signal. 

Included  within  the  analog  circuitry  is  an  amplifier  having  a  known  fixed 
(calibration)  resistor  in  its  feedback  loop.  The  signal  from  this  amplifier 
is  analyzed  simultaneously  with  the  rock  signal  and  provides  a  reference 
signal  of  known  amplitude  and  phase. 

Another  amplifier  is  provided  just  ahead  of  the  calibration  amplifier  to  pro¬ 
vide  a  means  of  easily  changing  the  gain  through  the  calibration  branch  with¬ 
out  changing  the  calibration  resistor  itself.* 
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The  raw  data  is  presently  extracted  as  printed  values  of  frequency,  series 
resistance,  series  reactance,  and  impedance  phase  angle.  The  method  of 
calculating  the  series  resistance  and  reactance  is  outlined  in  the  following. 
The  voltage  output,  eQ  (see  Figure  4),  from  the  amplifier  having  the  rock  as 
its  feedback  impedance,  is  given  by 


ZRei 


(1) 


Similarly,  the  output  e'  from  the  calibration  amplifier  is  given  by 


R'.R 
i  i 


e. 


(2) 


The  amplitude  of  the  complex  impedance  Z^  is  thus  given  by 


i  n 

KRi 


i 


(3) 


where  eQ  and  e'o  are  determined  from  a  Fourier  analysis  of  the  output 

voltages  e  and  e'  . 

o  o 

The  Fourier  analysis  also  provides  the  phase  angle  0  oetween  output  voltage 
and  current.  Thus  the  series  resistance  and  reactance  are  calculated  from 


Rs  =  |Zp  |  cos0,  and  Xg  =  JZp  |  sin  0 

A  least  square  circle  fit  is  also  automatically  made  for  the  R  and  X  data 

s  s 

and  displayed  along  with  the  R  and  X  data  on  a  CRT  screen.  In  addition, 

S  b 

all  data  are  permanently  stored  on  magnetic  tape  for  future  call  back. 
Detailed  description  of  the  computer  programs  and  electronic  circuitry  used 
in  this  work  can  be  found  in  reference  (1). 
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Figure  4.  Schematic  of  Experimental  Setup 
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B.  SAMPLE  PREPARATION 

With  the  exception  of  the  application  of  electrode  material,  all  measurements 
to  date  have  been  made  on  rocks  "as  received"  from  the  customer. 

The  electrodes  are  applied  to  the  flat  surfaces  of  the  cylindrical  samples. 

An  epoxy-based  conductive  adhesive  *  is  applied  in  a  thin  coat  and  allowed 
to  dry  for  at  least  two  days  at  room  temperature. 

Copper  leads  made  from  thin  sheet  stock  are  laid  on  top  of  the  electrode 
material  and  the  rock  is  clamped  tightly  between  two  plates  of  printed  circuit 
board  to  secure  even  electrical  contact  with  the  rock. 


$ 

Eccobond  Solder  V-91,  Emerson  &  Cumming,  Inc. 
Dielectric  Materials  Division,  Canton,  Massachusetts. 
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SECTION  III 

RESULTS  AND  DISCUSSION 

A.  IMPEDANCE  PARAMETERS  OF  DRESSER  BASALT 

A  disc  basaltic  sample  (I)  containing  Si02,  Fe^O^,  A^O^,  TiOg,  CaO  and 

MgO  in  the  percents  48.  42,  6.6,  15.23,  1.9,  8.  35  and  6.  4  percent,  respec- 

2 

tively,  and  measuring  0.  635  cm  in  length  and  6.  41  cm  in  cross  section  was 

supplied  to  Honeywell  by  the  Thermal  Fragmentation  Group,  Twin  Cities 

Mining  Research  Center,  U.S.  Bureau  of  Mines.  Silver  electrodes  were 

painted  to  the  circular  surfaces  of  the  rock  and  contact  to  the  electrodes  was 

made  with  copper  wires.  Replication  of  the  impedance  data  at  14  different 

frequencies  and  at  time  intervals  ranging  from  several  minutes  to  three  weeks, 

as  well  as  with  repeated  electrode  applications  to  the  cleaned  rock  surface, 

gave  fairly  reproducible  results.  The  standard  deviation  at  the  high  frequency 

of  2  KHz  was  3.  8  percent  of  the  mean  in  R  ,  and  2.  6  percent  of  the  mean  in 

s 

X  .  At  the  low  frequency  of  0.  05  Hz,  standard  deviations  of  3.  7  and  8.  1  of 
s 

the  mean  in  Rg  and  Xg,  respectively,  were  calculated. 

The  reactive  component  of  the  measured  impedance  is  plotted  in  Figure  5 
against  its  corresponding  resistive  component  at  a  series  of  frequencies.  The 
least- squares  circle  fit  to  the  data  was  determined  and  plotted  in  the  solid  arc 
of  Figure  5.  The  experimental  points  were  found  to  deviate  from  the  fitted 
arc  by  a  standard  deviation  of  2. 44  percent  of  the  fitted  radius.  Inspection  of 
the  data  in  Figure  5  suggests  a  value  of  16.  4  megohm  for  the  d-c  resistance, 
Rq,  of  this  basalt  sample.  A  rock  phase  angle,  0,of  29  degrees  was  also 
measured  for  this  sample. 

When  the  data  was  displayed  in  a  log  frequency  plot.  Figure  6  was  obtained. 
Curves  A  and  B  have  been  drawn  through  the  series  reactance  and  resistance 

Chemical  analysis  of  the  rock  samples  are  taken  from  a  recent 
report  by  Lindroth  and  Kranze  (reference  2). 
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Figure  5.  Argand  Circular-Arc  Plot  for  a  Basalt  Sample 
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Figure  6.  Dielectric  Dispersion  and  Relaxation  of 
Dresser  Basalt 
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data,  respectively.  Careful  examination  of  these  curves  indicates  three 
dielectric  relaxation  peaks.  The  main  peak  appears  at  a  frequency  of  21.  8 
Hz,  which  corresponds  to  a  turnover  frequency,  ^max,  of  136,  9  radians  per 
second,  and  to  a  relaxation  time,  t,  of  7.  3  milliseconds.  A  second  well- 
developed  relaxation  peak  appears  at  a  frequency  of  1.  7  Hz,  corresponding 
to  a  relaxation  time  of  93.  4  milliseconds.  A  third  relaxation  peak  appears 
to  be  anticipated  at  frequencies  less  than  0.  01  Hz.  In  agreement  with  these 
observations,  the  Rg  curve  shows  two  definite  dispersion  regions. 

The  first  dispersion,  appears  at  a  frequency  which  coincides  with  the 
first  and  main  relaxation  peak  in  the  series  reactance.  A  second  dispersion, 
c*2,  occurs  at  a  frequency  of  0.  05  Hz,  somewhat  displaced  from  the  second 
relaxation  peak.  A  third  dispersion  may  be  anticipated  at  frequencies  below 
0.  01  Hz.  In  conformity  with  the  adopted  notations  in  dielectric  dispersion, 
we  will  reserve  the  symbol  a  for  rock  dielectric  dispersion  in  the  low- 
frequency  range  (below  10  KHz).  The  symbol  0  will  be  used  for  dispersions 
in  the  medium  frequency  range;  i.  e. ,  between  10  KHz  and  100  megahertz. 
The  symbol  y  will  be  adopted  for  dispersions  at  frequencies  beyond  the  100 
megahertz.  Each  of  these  dispersions  may  split  into  more  than  one  sub¬ 
dispersion,  such  as  the  a^,  a^,  . . .  etc. ,  of  the  o-dispersion  investigated 
in  this  project.  Each  of  the  discovered  dispersions  should  correspond  to  a 
certain  relaxation  mechanism  of  a  substructure  "structon"  within  the  rock. 
Huggins  (Reference  3)  recommended  the  use  of  local  structural  groupings  as 
the  basic  structural  units.  He  used  the  term  "structon"  to  signify  a  specific 
type  of  atom(s),  with  specific  kinds  and  number,!  of  close  neighbors. 


B.  IMPEDANCE  PARAMETERS  OF  GRANITE 

A  sample  of  charcoal  gray  granite  from  St.  Cloud,  Minnesota,  was  supplied 
to  Honeywell  by  personnel  from  the  Thermal  Fragmentation  Group  of  the 
Twin  Cities  Mining  Research  Center.  The  sample  analyzed  (Reference  2) 
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G3.  5  percent  SiCg,  15.  6  percent  alumina,  4.  1  percent  CaO,  3.  6  percent  NagO, 
3.  6  percent  KgO,  4.2  percent  CaO,  2. 7  percent  FeO,  1.  8  percent  FegOg  and 
smaller  quantities  of  manganese  and  titanium. 


The  impedance  data  of  granite  followed  the  same  pattern  as  that  of  basalt, 

except  that  the  semicircular  arc  was  not  easy  to  close  at  the  low-frequency 

side.  Resistive  and  reactive  components  for  a  cylindrical  granite  sample 

2 

0.  148  cm  in  length  and  3.  67  cm  in  cross-sectional  area  described  a  segment 
of  a  circular  arc  when  plotted  in  an  Argand  diagram.  Computer  extrapolation 
of  this  arc  to  its  point  of  intersection  with  the  resistive  axis  yielded  a  value 

9 

of  4.  8  x  10  ohm  (about  5  gigaohm)  for  Rq.  The  preceding  basaltic  sample 
has  an  Rq  value  of  16.  4  megohm.  If  a  basalt  sample  of  the  same  dimensions 
of  granite  had  been  used,  its  Rq  value  would  be 


0.  635 

6.  41  j 

0.  148 

3.67) 

40.  5  megohm 


The  d-c  resistance  value,  Rq,  as  determined  by  extrapolation  of  the  impedance 
vector  to  zero  frequency,  is  therefore  about  a  hundred  times  higher  for  granite 
than  for  basalt  of  similar  shape  and  dimensions.  This  result  is  consistent 
with  the  fact  that  granite  contains  30  percent  more  of  the  covalently  coordinated, 
and  hence  nonconductive,  silica  tetrahedra.  Basalt  has  more  iron,  calcium 
and  titanium,  although  somewhat  less  sodium  and  potassium  than  granite. 


C.  IMPEDANCE  PARAMETERS  OF  QUARTZITE 

_3 

Sioux  quartzite,  with  a  bulk  density  of  2.  64  gram  cm  ,  was  used  in  the  fol¬ 
lowing  experiments.  The  source  location  of  this  type  of  quartzite  was  Jaspar, 
Minnesota.  Its  chemical  analysis  (Reference  2)  indicated  97.  84  percent  silica, 
0.  87  percent  AlgO^,  0.  81  percent  CaO,  0.2  5  percent  FeO  and  0.  2  7  percent 
^e2^3‘  6ran^e»  impedance  parameters  of  quartzite  were  very 

large,  in  the  gigaohm  range,  at  low  frequencies.  A  quartzite  disc  sample  of 
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length  0.  61  cm  and  cross-sectional  area  18.  8  cm  gave  an  extrapolated  R 

9  0 

value  of  4.  5  x  10  ohms.  When  corrected  to  the  same  dimensions  of  granite, 

this  value  should  be  modified  to 


4.  5  x  109 


0.  1 4  8  \ 
0.  61  j 


3.  67 
18.  8 


5.  6  x  109 


ohms 


Hence,  the  quartz  d-c  resistance  is  about  20  percent  higher  than  that  of 
granite,  and  still  more  than  two  orders-of-magnitude  higher  than  that  of 
basalt.  Thus  it  appears  that  1  ;  rock's  impedance  parameters  are  not  only 
affected  by  the  amount  of  silica  it  contains,  but  also  by  the  amounts  of  volatile 
alkali  oxides  within  the  rock. 


D.  EFFECT  OF  HUMIDITY  ON  ROCK  IMPEDANCE 

The  effect  of  moisture  content  has  been  examined  with  the  basaltic  sample  (I) 
and  the  data  are  shown  in  Figure  7.  Curve  A  in  this  figure  was  obtained  with 
the  rock  sample  as  received  from  the  customer  with  no  attempt  made  to  con¬ 
trol  the  rock  environment.  Curve  B  in  Figure  7  was  obtained  when  the  basalt 
sample  was  baked  with  electrodes  attached  in  a  110°C  oven  for  19  hours.  After 
baking,  the  sample  was  immediately  placed  in  an  airtight  desiccator  and 
allowed  to  cool  to  room  temperature  for  the  following  24  hours.  The  resulting 
measurements  indicated  a  significant  increase  in  impedance  at  frequencies  up 
to  500  Hz.  Above  500  Hz,  the  humidity  effect  on  impedance  parameters  is  not 
very  pronounced.  While  the  Rq  value  for  the  unbaked  basalt  was  16.  4  megohm, 
the  value  determined  when  the  rock  was  baked  is  about  45  megohm. 

Chemical  analysis  of  Dresser  basalt  (Reference  2)  indicated  the  presence  of 
0.  14  percent  water  in  the  rock  and  a  loss  on  ignition  (LOI)  of  2. 26  percent. 

By  contrast,  the  granite  sample  and  the  quartzite  sample  contained  0.  11  and 
0.  08  percent  water,  respectively,  and  showed  a  loss  on  ignition  of  0.  59  and 
0.  32  percent,  respectively.  Upon  ignition,  not  only  water  is  lost,  but  some 
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A  -  Rock  at  ambient  room  environment 
B  -  Dry  baked  rock  in  desiccator 


Figure  7.  Effect  of  Moisture  on  the  Impedance  Semicircle 
of  Basalt 
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of  the  voltatile  oxides  such  as  NagO  and  KgO  may  also  sublime.  When  these 
volatile  oxides  are  attached  to  the  nonvolatile  silica  tetral  edron,  they  are  in 
the  form  of  a  solid  solution  of  alkali  silicate,  and  the  alkali  will  thereby  be 
prevented  from  sublimation.  A  free  alkali  metal  oxide  unit  in  the  rock  will 
show  a  higher  LOI,  and  will  also  be  easily  leachable  with  water,  and  thus  will 
contribute  significantly  to  the  rock  conductivity.  Thus  the  alicali  content  of  a 
rock  as  shown  by  chemical  analysis  will  not  be  the  decisive  factor  in  deter¬ 
mining  the  rock  conductivity;  rather,  the  mode  by  which  the  alkali  oxide  enters 
the  rock  structure  will  play  the  predominant  role.  Because  of  its  higher  LOI, 
Dresser  basalt  appears  to  have  more  of  the  mobile  sodium  and  potassium  ions 
which  contribute  greatly  to  its  higher  conductivity. 


E.  EFFECT  OF  TEMPERATURE  ON  ROCK  IMPEDANCE 

Experiments  were  performed  on  basalt  to  determine  the  effect  of  temperature 
on  rock  impedance  as  a  preliminary  for  a  research  proposal  on  "Effect  of 
Frequency  and  Temperature  on  Rock  Dielectric  Parameters,"  which  was 
recently  submitted  to  ARPA  and  the  U.S.  Bureau  of  Mines.  The  basaltic 
sample  (II)  used  in  this  experiment  was  a  large  circular  disc  0.  54  cm  in 
thickness  and  5.  21  cm  in  diameter.  It  differed  from  the  Dresser  basalt  (I) 
used  previously  in  both  color  and  grain  size.  The  new  sample  was  darker 
in  its  greenish  tint  than  the  former.  The  change  in  texture  may  be  attributed 
to  different  phase  of  rock  formation  (Reference  4).  Impedance  has  been  mea¬ 
sured  at  both  room  temperature  (22°C)  and  the  melting  point  of  ice  (0°C).  At 
the  low  frequency  of  0.  1  Hz,  the  impedance  at  the  ice  point  is  considerably 
higher  than  that  at  room  temperature.  As  the  frequency  increases,  the 
temperature  effect  becomes  less  pronounced.  Table  1  gives  the  resistive  and 
reactive  impedance  components  at  room  temperature  and  at  nominal  frequen¬ 
cies  ranging  from  0.  1  to  2000  Hz.  The  variation  of  the  series  reactance  at 
room  temperature,  X  ,  with  log  frequency  is  shown  as  curve  A  of  Figure  8, 

b 

while  curve  B  shows  the  variation  of  the  rock  series  resistance,  R  .  Both 

s 


12288-IR1 


3-9 


Table  1.  Impedance  Data  on  Basalt  (II)  at  22°C 


LAST  RECORD  IS  T4-  2  4/25/71 

ZFEGO  MOD  I 


B6- 1  4/26/71 

ENTER  RCAL/UNI TS 

20. 

MEGS 


FRQ 

RES 

RAEC 

PHI 

0.09 

32.6651 

-3.5772 

-6.2  5 

0.50 

26.5789 

-4.5407 

-9.70 

0.99 

24.4720 

-5.3178 

-12.26 

2.00 

21.8083 

-6.1671 

-15.79 

5.02 

17.6153 

-7.0980 

-21.95 

7.95 

15.2352 

-7.2887 

-25.57 

9.96 

14. *585 

-7.2972 

-27.43 

14.97 

11.9454 

-7.1662 

-30,96 

20.02 

10.4838 

-6.9575 

-33.57 

30.03 

8.5419 

-6.5137 

-37.33 

40.16 

7.2568 

-6.1103 

-40.10 

50.10 

6.3452 

-5.7679 

-42.2  7 

99.50 

4.0188 

-4.5073 

-48.28 

200.2? 

2.3545 

-3.2393 

-53.99 

503.36 

1.1232 

-1.8585 

-58.86 

998.19 

0.6598 

-1.1528 

-60.22 

2022.06 

0.4074 

-0.6540 

-58.06 

17  SAMPLES 


CIRCLE  CENTER  =  16.9612  20.4344 

RADIUS  =  27.6828 

ERROR  0.88188  REDUCED  TO  0.58838 
IN  17  STEPS 
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I 


Figure  8.  Variation  of  Impedance  Parameters  with  Log 
Frequency  for  Dresser  Basalt  (II)  at  Room 
Temperature,  22°C 
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curves  indicate  a  turnover  frequency,  f  ,  of  9.  96  Hz,  which  corresponds 

max 

to  a  room  temperature  relaxation  time,  r,  of 

t  =  — 1 -  =  -  j- -  =  1 5.  9  milliseconds 

umax  max 


When  the  previous  experiment  was  repeated  at  the  melting  point  of  ice,  the 
data  in  Table  2  were  obtained.  Variation  of  the  impedance  parameters,  at 
room  temperature  with  log  frequency,  is  shown  in  Figure  9,  where  a  turn¬ 
over  frequency,  f*  ,  of  9.  96  Hz  is  determined.  At  0°C,  the  relaxation  time 
J  max 

is  given  by 


T 


/ 


1_ 
f 

max 


79.  5  milliseconds 


The  rock  d-c  resistance,  R  ,  at  0°C  is  191, 8  megohms,  which  should  be  com¬ 
pared  to  the  room  temperature  value  of  37.  3  megohm.  Thus  both  rock  d-c 
resistance  and  relaxation  time  decrease  significantly  with  rise  in  temperature. 
At  higher  frequencies,  the  temperature  effect  becomes  less  pronounced. 


According  to  Eyring's  theory  (Reference  6)  of  absolute  reaction  rates,  the 
turnover  or  characteristic  frequency  is  equal  to  the  universal  frequency 
(kT/h),  modified  by  a  free  energy  of  activation  term,  thus 


|  kTl  -AFt/RT 
max  '  h  I  e 


(4) 


where  k  is  Boltzmann’s  constant,  T  is  the  absolute  temperature,  AF*  is  the 
free  energy  of  activation  per  mole  of  the  relaxing  unit  within  the  rock  matrix, 
and  R  is  the  molar  gas  constant. 


If  it  is  assumed  that  the  silica  tetrahedral  unit  is  the  unit  structure  responsible 
for  the  relaxation  mechanism,  then  a  mole  of  the  relaxing  units  (Si04)  will  be 
92  grams  of  the  rock.  The  free  energy  of  activation  is  related  to  the  enthalpy 
of  activation,  AH,  per  mole  of  relaxing  units  by 
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Table  2.  Impedat.ee  Data  on  Basalt  (II)  at  0°C 


CRT  CALIBRATION 


DATA  TAPE  ON  UNIT  2  -  PUSH  START 
PAUS 

LAST  RECORD  IS  B6-1  4/26/71 
ZEEGO  MOD  I 


B7-1  4/26/71 

ENTER  RCAL/UNITS 

2D. 

MEGS 


FRQ 

RES 

RAEC 

PHI 

D.D9 

160.  3149 

-20.5281 

-9.40 

0.50 

11,0. 5509 

-36.2u07 

-18.16 

0.99 

92.3657 

-39.2870 

-23.04 

2.00 

72.3747 

-39.8985 

-28.87 

5.00 

47.1312 

-30.0387 

-37.41 

7.9b 

35.9375 

-32.2504 

-41.91 

9.95 

31.1345 

-30.1291 

-44.06 

20.02 

19.0096 

-22.9683 

-50.39 

30.08 

13.8184 

-18.9881 

-53.96 

40.23 

10.8011 

-16.3439 

-56.54 

50.05 

8.8994 

-14.4474 

-58.37 

99.50 

4.7105 

-9.4347 

-63.47 

200.27 

2.3633 

-5.7473 

-67.65 

504.20 

0.9996 

-2.7933 

-70.32 

998.19 

0.5700 

-1.5737 

-70.10 

2018.35 

0.3682 

-0.8424 

-66.40 

16  SAMPLES 


CIRCLE  CENTER  =  91.5041  87.5214 

RADIUS  =  129.8529 

ERROR  5.15521  REDUCED  TO  2.35518 
IN  11  STEPS 
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Figure  9.  Variation  of  Impedance  Parameters  with  Log  Frequency 
for  Dresser  Basalt  (II)  at  0°C 
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AF+  =  AHf  -  TAS*'  (5) 

where  AS*  is  the  entropy  of  activation.  Applying  Equation  (5)  into  (4),  one 
obtains 


f 


max 


kT\  A  S*/R  -AIlt/RT 
h  )  6 


(6) 


or 

’"<WT>  ■  '»*  +Tr,-Trt  (7> 

Normally  one  determines  fmax  at  a  series  of  temperatures  and  establishes  the 
validity  of  the  foregoing  equations  by  ascertaining  that  the  plot  of  In  (f  / T) 

against  (1/T)  is  linear.  From  the  slope  of  this  linear  plot  one  determines  the 
enthalpy  or  heat  of  activation,  A  H*.  thus 

Slope  =  (8) 

R 


and  from  the  intercept  of  the  linear  plot  with  the  ordinate,  one  can  determine 
the  entropy  of  activation,  A  S  ,  thus 

Intercept  -  In  {-^-J  +  (9) 

Applying  Equation  (6)  to  the  two  datum  obtained  in  this  preliminary  work,  and 
remembering  that  T  =  295°K  for  room  temperature  (22°C),  and  T  =  273°K 
for  the  ice  point,  then 


(f  /  T) 

max 

(f'  /  T) 

max 


or 


122  88- IH1 


3-15 


In 


-f  T/' 

max 

f_L 

r 

f/  T 

L  max  J 

R 

[t; 

‘  T 

(10) 


The  enthalpy  of  activation,  AH^, 


is  calculated  from  Equation  (10),  thus 


AH'  =  2.303R 


TT ' 
(T  -  TO 


log 


f  T ' 

max 

,  f'  T 

•-  max 


oao  i  no'?  295  x  273  , 

=  2.303  x  1.987  x - ^ -  log 

=  10,  700  calories  or  10.  7  Kcal. 


9.96  x  273 
2.0  x  295 


Hence,  it  is  estimated  that  the  activation  energy  for  the  relaxation  process  in 
basalt  is  about  11  Kcal  per  mole  of  relaxing  units,  or  about  0.5  ev  (1  ev  = 

23.  05  Kcal/mole) . 


Future  work  in  this  project  will  further  develop  this  temperature  effect  and 
utilize  it  to  compute  both  the  enthalpy  and  entropy  of  activation. 


F.  EFFECT  OF  ROCK  SIZE  AND  SHAPE  ON  ITS  IMPEDANCE  DIAGRAM 

Experiments  on  Dresser  basalt  samples  of  various  sizes  and  shapes  showed 
significant  variations,  particularly  at  the  low-frequency  end.  When  the  sample 
was  in  the  form  of  a  long  cylinder,  reasonable  agreement  was  obtained  in  the 
specific  impedance  values;  i.e.  ,  the  rock  specific  resistance  P  ,  and  specific 
reactance,  X.  These  values  were  obtained  from  the  measured  resistance,  R, 
and  reactance  ,  X,  by  multiplying  into  (A/d),  where  A  is  the  electrode  area 
and  d  is  the  length  of  the  sample,  or  the  distance  between  the  electrodes.  How¬ 
ever,  with  some  of  these  long  samples  the  impedance  Values  at  low  frequencies 
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were  so  high  that  measurements  were  difficult  at  lower  frequencies  (below 
100  Hz)  due  to  amplifier  saturation.  Upon  resorting  to  thinner  samples,  the 
electrode  polarization  effects  were  quite  substantial  and  constituted  a  signifi¬ 
cant  part  of  the  measured  impedance  especially  at  low  frequency.  Similar 
findings  were  reported  by  Takashima  and  Schwan  (Reference  4). 


G.  ELECTRODE  EFFECTS 

Thin  slices  were  cut  from  a  long  cylindrical  rock  sample.  The  same  electrode 
material  was  applied  to  both  sides  of  the  rock  sample  (with  the  longer  length 
d,  and  for  the  thinner  rock  sample  whose  length  was  d#).  The  total  impedance 
measured  with  either  rock  sample  is 

Z  =  Z  +  Z  (11) 

me 

where  Zm  is  the  measured  impedance,  Z  is  the  true  impedance  of  the  rock, 
and  Zg  is  the  electrode  impedance.  For  the  long  sample. 


and  for  the  thin  sample, 


where  ?  is  the  rock  impedivity  and  A  is  the  electrode  area.  Impedivity  § 
is  the  vector  summation  of  resistivity  p  and  specific  reactance  X  ,  thus 

§  =  P  +  jX  (14) 

Ry  separating  the  real  from  the  imaginary  variables,  identical  relationships 

to  those  derived  for  Z„  will  result  for  electrode  resistance  R  and  reactance 

e  e 

X  Substituting  from  (12)  and  (13)  into  (11),  one  obtains  for  the  long  sample 


12288-IR1 


3-17 


= 

m 

Mil  +Ze 

(15) 

and  for  the  thin 

sample 

Z  '  = 
m 

5  Ixl  +  ze 

(16) 

therefore. 


(17) 


where  a  is  the  ratio  of  the  lengths  of  the  two  rock  samples.  Solving  Equation 
(17)  for  ZQ,  one  obtains: 


<*Z'  -  Z 

m _ m 

(a  -  1) 


(18) 


Separating  Equation  (18)  into  its  real  and  imaginary  parts  gives 


»  R'  -  R 
m  m 

(«-  1) 


(19) 


and 


X 


* 

m _ m 

(a  -  1) 


(29) 


A  thin  slice  of  height  0. 152  cm  was  cut  from  a  basaltic  cylinderical  sample 

(III).  The  length  of  the  remaining  long  cylinder  was  4. 08  cm.  Hence,  a  =  26.  8. 

2 

Circular  silver  paste  electrodes,  each  of  area  3. 75  (cm  )  were  fastened  to 
each  side  of  the  two  cylinders.  The  impedance  data  measured  on  these  two 
rock  specimens  are  shown  in  columns  2,  3,  4  and  5  of  Table  3.  Re  and  Xg 
calculated  from  these  data  with  the  aid  of  Equations  (19)  and  (20)  are  given  in 
columns  6  and  7  of  Table  3. 
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By  subtracting  the  data  in  columns  6  and  7  point  by  point  from  the  corresponding 
measured  impedance  data,  and  mulitplying  each  resulting  number  into  (A/d), 
the  basaltic  specific  resistance  p  ,  and  reactance  X  ,  were  calculated  and 
recorded  in  columns  12  and  13  of  Table  3. 

A  check  on  the  valadity  of  this  technique  was  made  by  cutting  a  third  section  of 
the  basalt  sample  with  length  0,  150  cm.  Repeated  measurements  on  this  third 
sample  gave  data  which  agreed  closely  with  those  just  reported.  The  specific 
impedance  data  for  this  third  sample  were  calculated  and  found  to  be  within 
the  one  percent  error  from  those  recorded  in  Table  3.  This  agreement  sup¬ 
ports  our  technique  of  correcting  for  the  electrode  impedance  artifacts  and 
lends  credence  to  the  validity  of  our  impedance  measurements. 

From  here  on,  each  rock  sample  will  be  cut  into  two  unequal  lengths  and  the 
ratio  a  of  the  two  lengths  will  be  accurately  determined.  The  computer  pro¬ 
gram  will  be  modified  to  automatically  correct  for  electrodes  by  Equations 
(19)  and  (20)  before  the  final  data  of  p  and  x  will  be  displayed. 

In  a  similar  manner,  the  electrode  corrections  for  granite  were  applied  and 
the  data  in  Table  4  were  obtained.  The  plots  of  the  specific  reactance  V 
against  the  specific  resistance  P  at  a  series  of  frequencies  for  basalt,  granite 
and  quartzite  are  shown  in  Figure  10.  The  impedivity  of  each  of  these  three 
rocks  expressed  in  ohms  x  meter  is  thus  determined  at  frequencies  between  1 
and  2  000  Hz. 


H.  TRANSFORMATION  OF  IMPEDANCE  PARAMETERS  TO 
DIELECTRIC  PARAMETERS 

Ar  important  phase  of  this  research  is  to  provide  techniques  for  the  extraction 
oi  the  real  and  imaginary  parts  of  the  dielectric  constant  (€  )  from  their  imped¬ 
ance  counterparts,  namely  the  imaginary  and  the  real  parts  of  the  electric 
impedance,  Z. 
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p  (MEGOHM  -  METER) 


Figure  10.  Specific  Impedance  Parameters  for  A  -  Basalt; 
B  -  Granite;  and  C  -  Quartzite 
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i? 


t  , 


Thus, 


€  =  e'  -  Je 


// 


and 


Z  =  R  +  JX 

If  the  dielectric  is  symbolically  represented  by  a  simple  R-C  unit  (see  Appen¬ 
dix  A) ,  then  the  relation  between  X  and  R  on  the  one  hand,  and  e  "and  e'on 
the  other  hand  establishes  a  full  semicircle  in  the  Argand  diagram,  and  the 
Cole-Cole  plot  respectively.  The  circle  will  have  its  center  on  the  real  axis 
and  will  pass  by  the  point  of  origin,  i.  e.  ,  the  coordinates  of  the  circle  center 
will  be  (r,  0)  where  r  is  the  circle  radius.  Since  1;he  conformal  transformation 
(W  -  Z  *)  of  a  circle  results  in  a  circle,  a  plot  of  the  conductance  G  versus  the 
susceptance  uC,  instead  of  X  versus  R,  establishes  a  circular  arc. 

This  study,  has  developed  a  model  that  simulates  the  dielectric  dispersion  of 
a  real  system  (see  Appendix  B).  An  important  criterion  is  that  this  model 
possesses  a  single  time- constant  or  relaxation  time,  in  contrast  to  pre¬ 
vious  investigators  (References  7,  8,  9,  and  10)  who  used  a  distribution  of 
relaxation  times  to  account  for  the  dispersion  in  a  real  system  that  results  in 
a  circular  arc  with  a  depressed  center  in  the  Argand  or  the  Cole-Cole  diagram. 
In  the  following  procedure,  we  are  going  to  use  the  "bracketting"  technique 
which  relates  the  boundary  points  in  both  the  dielectric  and  the  impedance 
circular  arcs,  and  to  invoke  the  principle  that  the  same  relaxation  process, 
with  a  characteristic  time,  t  ,  is  responsible  for  both  the  dielectric  and  the 
impedance  arcs,  i.e. ,  >ve  assume  a  constant  relaxation  time. 

The  relaxation  time  for  a  dielectric  dispersion  process  following  the  Maxwell- 
Wagner  mechanism,  is  related  to  the  volume  resistivity  by  the  following 
equation  (Reference  11): 
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t  =  ee  p  (21) 

r  K 

_  12 

where  «r  is  the  permittivity  of  free  space  (8.  85  x  10  farad  per  meter)  and 
P  is  the  volume  resistivity  defined  by 

R  =  p  (^)  or  p  =  R  (j)  (22) 

where  A  is  the  rock  sample  cross-sectional  area  (electrode  area)  and  d  is  the 
thickness  of  the  rock  specimen  (the  interelectrode  distance). 

At  the  two  boundary  points  of  infinite  and  zero  frequency,  respectively,  the 
application  of  Equation  (21),  with  the  provision  that  the  value  of  C  decreases 
with  increasing  frequency,  leads  to 

t  =  t.trp0  .  erR0$> 

and 

T  =  €o  er  D®  -  cocrR-  <£> 

This  is  true  because  we  like  our  model  to  be  presented  with  a  single  relaxation 
time. 


(23) 


(24) 


The  dielectric  constants  at  infinite  and  zero  frequency  can  therefore  be  calcu¬ 
lated  from 


e 

r 


T 


and 


e 

r 


€ 

o 


(25) 


(26) 
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The  relaxation  time  can  be  related  to  our  model  parameters  (see  Appendix  B), 
in  the  following  manner: 


t  =  (Rq  -  Rg,)  Cp  Vl  +  cotan*'  0 


(27) 


where  0  is  the  phase  angle  derived  from  the  impedance  diagram  as  follows: 


=  sin  _1  [  Ro  ~  R» 

L  2r 

where  r  is  the  radius  of  the  circle  whose  arc  represents  the 
The  Cole-Cole  parameter,  a,  is  related  to  the  phase  angle. 


(2  8) 

Argand  diagram. 
0  ,  as  follows: 


0  =  J  (1  -  a)  (29) 

In  the  ideal  case  of  a  Debye  dielectric,  0  = and  a  ~  In  general,  a  is  a 
fraction,  usually  close  to  0.  5. 


The  capacitance,  Cp  needed  in  Equation  27  for  calculating  the  relaxation 
time  can  be  obtained  from  the  impedance  data  as  follows: 

-X 

cp  =  _ _  (30) 

“(Ru  -  R->2  +  xj 

Plotting*  Cp  as  a  function  of  frequency,  a  sharp  decrease  occurs  between  0.01 
and  100  Hz,  which  is  followed  by  a  plateau  between  100  Hz  to  2,000  Hz.  The 
average  value  of  this  plateau  is  considered  to  represent  Cp. 

Using  the  phenomenological  equations  of  the  complex  dielectric  constant 
(Referenes  7,  12). 
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e  =  e 


«:  -  e„ 
o  ® 


®  i  -a 

1  +  (Jut)1 


=  e'  -  Je" 


(31) 


Equation  (31)  is  separable  into  its  real  and  imaginary  parts  using  the  identity: 

(32) 


.a 

3  =  e 


in.  a 
7 


to  give 


e  -  efl 


£'  =  £„+  -2— 
CO  2 


1  - 


sinh  (l-o)  S 


cosh  (1 -o)  S  +  cos  (-rr-) 


T'J 


(33) 


and 

■9  ~  %a  ^  COS  (ir) 

e //  _  2  o  t»° _ 2 

cosh  (1  -o')  S  +  sin  (^jp) 
where  S  is  given  by 

S  =  loge  (ut) 

The  plot  of  e"  (w)  as  a  function  of  e '  (u)  as  calculated  from  Equations  (33)  and 
(34)  gives  the  Cole-Cole  diagram. 

The  necessary  operational  steps  in  accordance  with  the  foregoing  equations 
were,  therefore,  added  to  the  computer  algorithm  to  extract  the  dielectric 
permittivity  data  from  the  specific  impedance  parameters.  Provision  for 
correcting  the  measured  impedance  parameters  for  electrode  effects  were  also 
included  in  the  program.  The  following  subsection  summarizes  the  dielectric 
data  obtained  with  this  procedure. 


(34) 


(35) 
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I.  COLE -COLE  DIAGRAMS  -  DIELECTRIC  AND  IMPEDANCE 
PARAMETERS  OF  DRESSER  BASALT 

Basaltic  sample  (ill)  was  obtained  in  the  form  of  a  cylinder  with  diameter  2.17 
cm.  Thin  slices  that  measured  0.30  and  0.61  cm, in  length  were  cut  from  this 

A 

sample.  The  ratio  (^-)  for  the  thinner  slice  was  0.  1232  meter  and  for  the 
thicker  slice  was  0.0618  meter.  The  electrode  impedances  were  calculated 
with  the  aid  of  Equations  (19)  and  (20)  for  the  two  slices.  The  electrode 
impedances  were  subtracted  point  by  point  from  the  corresponding  measured 
impedance  values.  The  resulting  data  were  then  multiplied  into  (A/d)  to  give 
the  specific  resistance  and  reactance  of  basalt.  Photostatic  copies  of  the 
computer  output  for  the  two  basaltic  slices  are  shown  in  Tables  5  and  (i. 

Dielectric  permittivity  at  each  frequency  was  calculated  from  the  model 
parameters  as  described  in  subsection  H.  The  Cole-Cole  plots  resulting  from 
these  data  are  shown  in  Figures  11  and  13  for  the  0.30  and  0.61  cm  slices  of 
basalt,  respectively.  These  figures  can  be  considered  to  be  in  satisfactory 
agreement  within  experimental  error.  The  dielectric  dispersion  curves  for 
basalt  are  shown  in  Figures  12  and  14  for  the  0.30  and  0.61  cm  slices  of 
basalt.  The  relaxation  times  calculated  from  either  the  maximum  in  dielectric 
loss,  e",  or  the  inflection  point  in  the  e'  curve  were  found  to  be  16. 6  and  13.4 
milliseconds  for  the  two  basalt  slices.  Since  the  experimental  measurements 
on  both  slices  were  completely  independent  of  each  other,  the  derived  values 
of  relaxation  time  are  considered  to  be  in  reasonable  agreement. 

Table  7  summarizes  the  impedevity  and  dielectric  permittivity  data  on  the 

two  basalt  slices.  The  dc  resistivity  of  basalt,  P  ,  is  estimated  as  1.27  x 
8  8  o 

10  and  1.00  x  10  ohm  meter  for  the  thin  and  thick  slices,  respectively. 

These  resistivity  data  appear  to  be  internally  consistent  with  each  other  and 

8 

are  in  agreement  with  the  value  1.26  x  10  ohm  meter  reported  in  Parkhomenko 
(Reference  13)  for  dry  basalt.  The  agreement  is  surprising  in  view  of  the 
anticipated  differences  between  our  Dresser  basalt  and  Parkhomenko's  most 
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Table  5.  Impedance  and  Dielectric  Data  on  Basalt  (III) 
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-41.8617 

412.9874 

0.601 

227.0444 

-224.0116 
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RS,XS 
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RESULTS-  CIRCLE  CENTER  s 

511.1931  370.7580 

RADIUS  = 

638.4723 

FIT  = 

3.4752 

R0= 

1030.9866 

P1NF-  -8.6005 

MODEL 

PARAMETERS 

R 1  =  1545.4873/FRO 

RP  =  1039.5872 

CP  *  0. 1701E-03 


CP=  0.130032E-04 


TAU-  16 .O0434MSEC 


E  0  = 1  7  7  1.  1357 


E1NF=  14.77477 


FRQ 

EP 

EPP 

0.099 

-0.168166E 

04 

-0.851248E 

02 

0.150 

-0, 16  5  7  8  IE 

04 

-0, 106640E 

03 

0.201 

-0 . 163703E 

04 

-0. 12474 2E 

03 

0.403 

-0. 157142F 

04 

-0. 178948E 

03 

0.601 

-0. 152152F 

04 

-0.217142E 

03 

0.801 

-0.  14794  IE 

04 

-0.247313E 

03 

1.002 

-0 • 1442S7E 

04 

-0.272144E 

03 

2.012 

-0.13032  IF 

04 

-0.352604F 

03 

4.022 

-0. 112952E 

04 

-0,42212  IE 

03 

0.009 

-0.101618E 

04 

-0.448484E 

03 

10.052 

-0 . 8640070 

03 

-0.459785E 

03 

20.072 

-0.6O2553E 

03 

-0 . 4  32068F 

03 

60.096 

-0.393816F 

03 

-0.320979K 

03 

100.004 

-0.798 12 5E 

03 

-0.2O1763E 

03 

20  1.  362 

-0. 199427C 

03 

-0. 190306E 

03 

003.440 

-  0 . 994948E 

02 

-0.107468E 

03 

1009.174 

-0,O9742:r 

02 

-0.807935E 

02 

1509.434 

-0. 518329E 

02 

-0.642934E 

02 

2025 , 783 

-0,411 46 7F 

02 

-0.542893F 

02 

12288-IR1 


3-27 


Table  6.  Impedance  and  Dielectric  Data  on  Basalt  (III). 
0.  61  cm  slice. 
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Figure  11.  Cole- Cole  Plot  of  Dresser  Basalt.  Data  on 
0.  61  cm  slice. 
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Figure  12.  Dielectric  Permittivity  of  Dresser  Basalt  as  a 
Function  of  Frequency  (a  -  Real  Part  of  Dielec¬ 
tric  Constant,  o';  b  -  Imaginary  Part,  e'/). 
Data  on  the  0.  61  cm  slice;  relaxation  time, 
t  -  13.4  millisecond. 
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Figure  14.  Dielectric  Permittivity  of  Dresser  Basalt  as  a 
Function  of  Frequency  (a  -  Real  Part  of  Dielec¬ 
tric  Constant,  e  ;  b  -  Imaginary  Part,  z"). 
Data  on  the  0.  30  cm  slice;  relaxation  time, 
t  =  16.  6  milliseconds. 
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Table  7.  Electric  and  Dielectric  Parameters  of 
Dresser  Basalt 


Parameter 

0.  30  cm  Slice 

0.61  cm  Slice 

dc  resistivity,  pQ.  ohm  meter 

1.27  x  108 

1.00  x  108 

infinite  frequency  resistivity, 

p  ,  ohm  meter 

00 

1.06  x  106 

1.01  x  106 

model  parallel  resistance,  R^. 
ohm 

8.  89  x  1014 

-5 

6.  87  x  1014 

-6 

model  capacitance,  C  ,  farad 

1. 3  x  10  ° 

6. 8  x  10 

relaxation  time,  t  ,  second 

16.6  x  10"3 

13.4  x  10"3 

zero  frequency  dielectric 
constant,  eQ 

1.77  x  103  . 

1.50  x 103 

infinite  frequency  dielectric 
constant,  e 

00 

14.  8 

14.  8 

maximum  dielectric  loss, 

Gmax 

4.6  x  102 

4.2  x  102 

poG~  V  second 

16.4  x  10“3 

13. 1  x  10_3 

pcoGo  er  second 

16.6  x  10'3 

_3 

13.2  xlO 

probably  Russian  basalt  (origin  not  given).  Table  7  also  shows  that  the  zero 

Q  Q 

frequency  dielectric  constant  of  basalt  is  1.77  x  10  and  1.50  x  10  for  the  thin 
and  thick  slices.  The  internal  agreement  in  this  case  is  also  acceptable. 
Parkhomenko  (Reference  14)  reports  that  at  low  frequencies,  10^  to  10^ 

o  c 

Hz,  the  dielectric  constant  assumes  very  large  values  (10  to  10  ). 
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Keller  (Reference  15)  observed  that  the  product  of  the  dielectric  constant  at 
low  frequency,  eQ,  into  the  resistivity  at  low  frequency;  p  Q  ,  is  nearly  a 
constant  characteristic  of  a  particular  type  of  rock.  For  rhyolite  and  basalt, 
the  average  value  of  pQ  is  0.G3  sec  with  the  range  of  19  values  from  0.25 
2G  sec.  Average  values  of  this  product  varied  from  1.  4  x  10"^  for  basic 
igneous  rocks  such  as  gabbro  and  chromite  to  10.  8  for  the  hematite  ore  from 
Minnesota,  which  contains  appreciable  quantities  of  electronically  conducting 
materials.  The  data  in  Table  7  also  show  that  the  products  p  e*  e  and 
P<DC0er  are  equal  to  the  relaxation  time  for  both  the  thin  and  thick  samples  of 
basalt. 

Similar  experiments  to  those  described  above  were  run  on  a  cylindrical 

A 

quartzite  sample  of  2.22  cm  diameter,  0.  955cm  in  length  and  an  (-g)  ratio 
of  0.  0496  meter.  The  computer  output  data  on  this  sample  are  shown  in 
Table  8.  From  the  present  data,  one  can  calculate  the  d-c  resistivity  of 
quartzite,  pQ,  from  the  observed  value  Rq  of  31737.  0  Megohm,  thus 

A.  10  q 

PD  =  Rq  ("j)  =  3.  17  x  0.041  x  10  =  1.29  x  10  ohm-meter 

The  relaxation  time  of  this  quartzite  sample  is  also  calculated  to  be  GO.  3 
milliseconds.  The  real  component  of  the  dielectric  constant  at  infinite  frequency, 
,  is  determined  as  5.29.  Also,  the  real  component  of  the  quartzite  dielectric 
constant  at  the  limit  of  zero  frequency,  ,  is  found  to  be  3.  95  x  10  . 

Further  trials  are  currently  in  progress  to  improve  the  model  used  to  derive 
the  rock  electrical  parameters  and  the  steps  by  which  impedance  data  are 
transformed  to  dielectric  data. 
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Table  8.  Impedance  and  Dielectric  Data  on  Quartzite 
0.  30  cm  slice. 
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SECTION  IV 

TECHNICAL  REPORT  SUMMARY 


The  research  program  described  in  this  report  was  initiated  six  months  ago 
at  Honeywell  Systems  &  Research  Division  with  the  objective  of  applying 
dielectric  relaxation  theories  to  rocks,  and  to  utilize  the  dielectric  parameters 
and  the  closely-related  electrical  impedance  parameters  to  obtain  significant 
correlates  with  rock  geophysical  characteristics  and  ultimate  structure.  The 
information  gathered  will  be  pertinent  to  rock  fragmentation,  excavation  and 
underground  tunnelling  problems.  The  degree  of  attenuation  of  magnetic 
fields  from  a  flat-lying  loop  in  an  underground  tunnel  depends  on  the  electri¬ 
cal  resitivity  of  the  overlying  rock  and  hence  electromagnetic  communication 
with  underground  tunnels  and  mines  will  be  greatly  improved  from  an  accu¬ 
rate  knowledge  of  rock  impedance  and  dielectric  properties.  The  data 
gathered  will  also  be  relevant  to  the  presence  of  underground  water,  entrapped 
water  and  gases  in  rocks,  and  similar  factors  needed  ahead  of  excavation 
and  underground  tunnelling.  To  achieve  these  goals,  accurate  and  rapid 
methods  of  displaying  rock  electrical  impedance  (both  resistive  and  reactive) 
and  dielectric  permittivity,  will  first  be  developed. 

Our  present  novel  techniques  of  impedance  data  display  permits  direct  evalu¬ 
ation  of  the  dc  resistance,  Rq  (which  is  directly  related  to  dielectric  permit¬ 
tivity  at  zero  frequency,  eQ);  the  resistance  at  infinite  frequency  (which  is 
directly  related  to  the  permittivity  s®);  and  the  impedance  (and  complex 
dielectric  constant)  at  any  frequency  in  between.  Hence,  as  a  first  step, 
our  existing  impedance  equipment  with  its  on-line  computer  capabilities  was 
modified  to  measure  the  complex  impedance  of  rocks.  A  brief  description 
of  the  new  equipment  is  given  in  Section  II  of  this  report.  It  is  highly 
desirable  to  extend  the  capability  of  this  equipment  (beyond  the  kilohertz 
frequency  range).  Extension  to  the  mega-  and  gigahertz  range  has  been 
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recommended  in  Honeywell  proposal  1D-E-3,  "Effect  of  Frequency  and 
Temperature  on  Rock  Dielectric  Parameters",  which  was  recently  submitted 
to  the  contracting  agency. 

The  complex  impedance  of  rock  samples  at  a  range  of  frequencies  that 
brackets  their  turnover  frequency  would  yield  valuable  information  on  the 
relaxation  time(s)  of  the  silica  tetrahedra  and  other  structural  groups  within 
the  rock,  and  the  presence  or  absence  of  conductive  minerals.  In  general, 
a  Cole -Cole  (Reference  7)  circular -arc  plot  would  be  obtained  in  the  complex 
dielectric  diagram.  With  complex  structures,  the  resulting  figure  may  be 
analyzed  into  a  eries  of  nearby  semicircular  arcs,  each  describing  a  given 
relaxation  mechanism  with  some  interaction  coefficients  among  the  various 
aggregates.  The  same  Cole-Cole  plot  can  be  obtained  in  a  complex 
impedance  diagram,  sometimes  called  Argand  diagram.  Sinbel  (Khalafalla, 
Reference  16)  showed  that  the  plot  of  the  series  reactance,  Xs,  against  the 
series  resistance,  R  (both  determined  at  a  given  frequency),  in  a  complex 
system  should  yeild  a  semicircular  arc.  The  intersection  of  this  circular 
arc  with  the  r  ial  (resistive)  axis  will  define  the  dc  resistance,  Rq  (at  the 
farthest  right  end).  Any  point  on  the  arc  will  define  the  impedance  radius 
vector  with  both  its  reactive  and  resistive  components.  The  results  shown 
in  Figures  1,  2,  3  and  5  in  Section  I,  and  7,  8,  11  and  13  in  Section  III  amply 
support  the  validity  of  our  concepts.  Impedance  data  were  gathered  on 
basalt,  granite  and  quartzite,  and  despite  the  wide  variation  in  electrical 
properties  of  these  rocks,  we  were  able  to  fit  their  data  to  appropriate 
circular  arcs. 

Relaxation-times  and  rock  characteristic  frequencies  can  be  easily  deter¬ 
mined  by  displaying  the  X  data  and  R  data  as  functions  of  frequency  in 
semilogarithmic  graphs.  Theory  indicates  that  X  (and  also  the  imaginary 
part  of  dielectric  permittivity,  e",  or  the  dielectric  loss)  exhibit  a  maxi¬ 
mum  value  at  the  turnover  or  characteristic  frequency.  On  the  other  hand, 
the  resistive  component  of  the  impedance,  R  (and  also  the  real  part,  e' ,  of 
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the  complex  dielectric  permittivity)  will  undergo  a  sudden  dispersion  at  the 
turnover  frequency.  Our  experimental  data  shown  in  Figures  6,  8,  9,  12 
and  14  in  Section  III  are  in  full  agreement  with  these  predictions.  The  appli¬ 
cation  of  these  fundamental  impedance  and  dielectric  theories  to  rocks  is  rather 
novel,  and  is  only  possible  because  of  the  accuracy  and  rapidity  of  our  technique 
of  data  acquisition. 

Electrode  polarization  effects  are  always  encountered  in  electrical  measure¬ 
ments  of  complex  systems.  A  method  is  developed  in  Section  III,  subsection  C 
to  correct  for  electrode  polarization  by  cutting  thinner  slices  of  the  same  rock. 
The  method  was  applied  with  success  on  the  impedance  data  of  basalt  and 
granite.  The  electrode  effects  are  significant  for  thin  rock  samples  where 
they  may  constitute  a  substantial  segment  of  the  measured  impedance.  While 
the  electrode  corrections  may  be  neglected  for  longer  rock  samples,  sometimes 
impedance  measurement  may  be  exceptionally  difficult  because  of  amplifier 
saturation  at  the  very  high  impedance  levels.  It  is  concluded  that,  whenever 
possible,  electrode  polarization  effects  should  be  determined,  and  corrected 
for,  by  repeating  the  impedance  measurement  on  a  thin  slice  of  the  same 
rock  sample. 

The  effect  of  temperature  and  pressure  on  relaxation  time  and  characteristic 
frequency  can  be  used  to  calculate  the  activation  energy,  activation  entropy, 
and  activation  volume  of  the  aggregate  interactions.  This  new  set  of  data 
can  yield  further  information  on  the  ultimate  structure  and  petrogenesis  of 
rocks. 

These  new  areas  of  endeavor  constitute  the  Honeywell  long-range  research 
goals  and  were  submitted  in  a  new  proposal  1D-E-3  to  the  contracting  agency. 
To  validate  the  concepts  presented  in  this  proposal,  we  performed  some  pre¬ 
liminary  measurements  on  basalt  at  two  temperatures.  The  results,  given  in 
Section  HI,  subsection  G,  of  this  report,  indicate  a  significant  temperature 
effect  on  both  the  relaxation  time  and  the  dc  resistance  of  basalt.  An 
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activation  energy  of  about  0.  5  ev  for  the  relaxation  of  structural  units  within 
basalt  was  estimated.  Further  work  on  this  interesting  temperature  effect 
will  hopefully  be  conducted  in  the  second  phase  of  the  project;  i.  e. ,  in  the 
following  year. 

A  new  technique  was  developed  for  the  extraction  of  the  dielectric  parameters 
from  the  measured  impedance  parameters  after  electrode  correction  and 
size  normalization  was  adopted.  The  method  depends  on  "bracketting"  the 
circular  arc  in  both  the  Argand  and  the  Cole -Cole  plots  and  on  the  assumption 
that  the  single  relaxation  time  must  be  the  same  for  both.  The  computer 
algorithm  has  been  extended  to  derive  the  dielectric  loss,  eM,  and  the  real 
part  of  dielectric  constant,  s',  at  various  frequencies  from  the  measured 
specific  impedance  parameters.  The  results  on  basalt  agree  reasonably  well 
with  those  reported  in  literature  and  are  described  in  Subsections  H  and  I  of 
Section  III. 

Research  efforts  in  the  following  six  months  will  be  devoted  to  developing  and 
improving  the  electrical  models  that  are  capable  of  simulating  the  complex 
impedance  of  rocks,  and  to  apply  these  models  to  the  gathered  data.  Further 
emphasis  will  be  placed  on  investigating  the  effects  of  rock  size,  rock  shape 
and  dimensions,  as  well  as  degree  of  dryness  and  humidity  on  the  generated 
impedance  data. 
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APPENDIX  A 

ELECTRICAL  MODELS  TO  SIMILATE  ROCK  IMPEDANCE 


Impedance  data  at  various  frequencies  usually  describe  a  circular-arc  plot 
when  displayed  in  a  cartesian  diagram  with  the  series  reactance  as  ordinate 
and  the  series  resistance  as  abscissa.  The  display  is  never  a  full  semi¬ 
circle,  but  is  a  part  of  a  circle  whose  center  does  not  lie  on  the  real  axis 
(resistance),  and  is  substantially  below  it.  The  phase  angle  is  defined  as  half 
of  the  angle,  subtended  by  the  circular-arc  at  the  center  of  the  circle,  i.  e, , 
between  the  two  radii  of  the  circle  defining  the  two  points  of  intersection  of 
the  arc  with  the  resistive  axis.  An  ideal  dielectric  dispersion  model  system 
is  easily  simulated  by  a  parallel  R-C  unit,  and  can  be  shown  to  give  a  full 
semi-circular  plot  upon  transformation  to  an  iso -impedic  series  R-C  unit, 

A  major  objective  of  this  research  is  to  seek  for  correlations  between  one  or 
more  of  the  rock  impedance  parameters  and  its  geophysical,  structural,  as 
well  as  genetic  characteristics.  If  significant  correlations  are  discovered, 
then  the  former  parameters  can  be  taken  as  a  quantitative  predictor  of  the 
latter  measures.  As  a  first  step  to  achieve  this  goal,  an  adequate  model 
capable  of  simulating  the  rock  impedance  data  should  be  developed.  This 
constitutes  an  important  topic  of  the  proposed  investigation. 

Sinbel  (Khalafalla)  (Reference  16)  derived  various  analytical  proofs  of  the 
semicircular  arc  in  an  attempt  to  describe  the  relationship  of  the  equivalent 
series  reactance  to  the  corresponding  resistance.  The  derivations  were 
based  on  a  parallel-to-series  transformation  of  electrical  models  with  a 
single  time -constant.  In  the  simplest  case,  it  was  shown  that  the  trans¬ 
formation  from  a  parallel  R-C  unit  (with  frequency  independent  components) 
to  an  adjustable  series  R-C  unit  (with  the  same  total  impedance)  results  in  a 
full  semicircular-arc  plot  in  the  Argand  diagram,  "’"he  locus  of  the  series 
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reactance,  X  ,  when  plotted  against  the  series  resistance,  R  ,  follows  the 
analytical  equation  of  a  circle  of  radius,  y  Rp,  where  Rp  is  the  parallel 
resistance  assumed  to  be  a  constant  quantity  characteristic  of  the  system,  thus 


X  2  +  [R  -  4-  R  ]2 

s  s  2  p 


1  R  2 

T  RP 


(Al) 


The  center  of  the  semicircular  plot  should  then  have  the  coordinates  (0,  -g-  Rp) 
and  hence  is  located  on  the  real  or  resistive  axis.  This  seems  to  represent 
a  very  idealized  situation.  In  real  systems  the  semicircular  arc  observed 
experimentally  is  usually  translated  vertically  downwards  so  that  its  center 
has  the  coordinates  a  and  -b,  and  accordingly  follows  an  equation  of  the  form 


[X  +  b]2  +[R_  -  aJ2  =* 


(A2) 


where  a,  b  and  c  are  constants  related  to  X  and  C  , 

r  r 


The  assumption  of  a  leaky  condenser  was  made  by  connecting  a  small  resis¬ 
tance,  r,  to  the  parallel  capacitor,  Cp,  in  the  idealized  parallel  R-C  model. 

This  resulted  in  relationships  between  X  and  R  of  the  form 

s  s 


Rs'r 


R 


(r  + 


i  P 


R_ 


4(r  +  Rp)2 


(A3) 


This  again  represents  an  equation  of  a  circle  whose  center  lies  on  the  real 


axis  at  a  distance  of  R  (r  +4-  R  )/(r  +  R  )  from  the  origin  and  whose  radius 
l  P  c  P  P 

is  given  by-^-  R  /(r  +  R  ), 
c  P  P 


Several  other  alternatives  were  tried  (Reference  16)  by  attaching  resistors 

and  capacitors  to  the  fundamental  R-C  unit  at  various  places.  However, 

P  P 

the  ensuing  circular  arc  plot  did  not  in  any  case  possess  the  limiting  and 
boundary  values  expected  from  the  electrical  model.  The  developed 
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mathematical  equations  were  unable  to  show,  for  example,  that  the  circular 
arc  plot  would  intercept  the  real  axis  at  points  corresponding  to  r  and  (r  +  It  ) 
as  would  be  expected  at  the  limits  of  infinite  and  zero  frequencies,  respectively. 
The  futility  was  evident  of  describing  a  complex  system  by  an  electrical  analog 
with  conventional  frequency  independent  components  and  a  single  time -constant. 

A  distribution  of  relaxation  times  can  in  principle  yield  expressions  Lof  which 
Equation  (31)  in  Section  II  is  a  special  case]  that  are  capable  of  explaining  the 
observed  results.  This  distribution  is  simulated  by  an  infinite  series  of 
parallel  R-C  units  with  various  time -constants.  Unlike  electrical  models 
with  a  single  time-constant,  these  mathematical  simulations  will  probably 
fail  to  provide  a  one-to-one  correspondence  between  the  rock  impedance 
parameters  and  its  geophysical  characteristics. 
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APPENDIX  B 

ELECTIUCAL  ANALOG  WITH  ONE  FREQUENCY 
DEPENDENT  RESISTANCE 


An  electrical  model.  Figure  HI  in  which  a  frequency  dependent  resistance, 
rj,  is  shunted  across  the  condenser,  C^,  is  capable  of  describing  a  circular 
arc  in  the  series  domain  with  a  vertically  displaced  center  (Figure  B2),  A 
condition  to  be  imposed  on  the  resistor,  y,  is  that  its  value  changes  inversely 
with  the  frequency  such  that  y  =y-,  where  g  is  a  constant  and  f  is  the  frequency. 
The  total  impedance  of  this  model  is  given  by 


R 


jr,X 
J  lp 

p  ri  +>\ 
>r\\ 

uP  +  ri  +  )xp 


+  r.. 


riRpxp(ri+V 

2-2  ,  „  ,2 

rinp  +  xp<ri +  V 


+  1 


rl  RnXn 
1  p  p 


r2R2  +  X2  (r  +  R  )2 
1  p  pi  p 


+  r 


2(B1) 


=  Rs  +  jXs* 


From  which  one  obtains 


R 


r2  + 


^y'p(ri  +  v 

qup  +  xp  (ri  +  R/ 


(B2) 
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Figure  Bl.  Electrical  Model  That  Produces  a  Circular 
Arc  Plot 


Figure  B2,  Circular  Arc  Plot  Between  Equivalent  X 

and  R  of  the  Model  s 

& 
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ncl 


X 


•>  '> 
i'T  ir  x 
j> . P 


•>  2  •'  2 

mil  +  X  (i-  +  R  ) 

!  p  pi  p' 


(B3) 


nking  the  following  representative  values  for  the  electrical  parameters, 


K  =  r0  =  300  ohm 

p  2 


C  =1  microfarad 
P 

g  =  150,000  ohm  per  second,  such  that  r^  =  “  300  ohm  at  a 

frequency  of  500  Hz,  one  is  able  to  calculate  values  for 

R  and  X  at  various  frequencies.  These  calculations 
s  s 

are  given  in  Table  Bl. 


Table  B 1 .  Equivalent  Series  Resistance  and  Reactance 
for  Electrical  Model  of  Figure  Bl 


f 

11  z 

rl 

ohm 

XP 

ohm 

Rs 

ohm 

X 

s 

ohm 

1 

150,  000 

159,  236 

599 

0.  6 

10 

15,  000 

15,  924 

594 

5.4 

100 

1,  500 

1,  592 

544 

38.  3 

3  00 

500 

531 

467 

58.  9 

500 

300 

318 

423 

57.  8 

700 

214 

227 

396 

52.  8 

1.  000 

150 

159 

372 

4  5.  1 

10,  000 

15 

16 

308 

7.  1 
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The  plot  of  X  against  R  as  shown  in  Figure  B2  is  found  to  describe  a 
s  s 

circular-arc  with  a  depressed  center  and  whose  boundary  values  satisfy  the 

limiting  values  required  by  the  electrical  models  of  the  circuit  in  Figure  Bl. 

Thus,  at  infinite  frequency,  both  r.  and  X  become  zero,  and  the  total 

impedance  of  the  network  equals  r i.e.,  300  ohm.  The  first  intersection 

point  of  the  circular-arc  with  the  Rg  axis  is  also  found  to  be  300  ohm.  At 

zero  frequency,  the  limiting  impedance  of  the  network  becomes  (R  +  r„) 

p  z 

=  600  ohm. 

It  appears,  therefore,  that  the  electrical  model  of  Figure  Bl  gives  a  realistic 
analog  to  a  system  displaying  a  circular-arc  plot  in  the  Argand  diagram.  The 
limits  of  the  rock  impedance  at  zero  frequency  would  be  given  by  (r9  +  R  ) 

<2  p 

and,  at  infinite  frequency,  by  A  direct  estimation  of  r 2  and  Rp  is  thus 
possible  by  extrapolating  the  circular-arc  to  intercet  the  real  axis.  The 
unusual  resistor  in  this  model  is  r1  =  where  g  is  a  constant  given  by 

g  =  2nKC  <B4> 

r 


where  K  is  a  constant  related  to  the  phase-angle,  and  C  is  the  frequency 
■  P 

independent  rock  capacitance  which  can  be  taken  as  a  constant. 


The  rock  phase-angle,  0,  would  only  be  dependent  on  the  values  of  Cp  and 
rj  and,  since  these  two  components  are  in  parallel,  then 


0 


tan 


tan 


(B5) 
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and,  hence, 

K  =  cotan  0,  or  0  =  cotan"*K. 


(B6) 


The  angle,  0,  whould  also  be  identical  with  the  experimentally  measured 
phase-angle  between  the  vertical  through  the  semicircular -arc  center  and  the 
line  joining  it  to  either  of  its  intercepts  with  the  real  axis.  Thus,  experimentally 
determining  the  rock  phase-angle,  0,  in  each  individual  case  would  enable  a 
determination  of  the  constant  K  required  to  calculate  r^. 

To  calculate  the  turnover -frequency,  one  can  make  use  of  the  fact  that  it  is 

the  characteristic  frequency  which  maximizes  X  .  Upon  differentiating 

s 

Equation  (B3)  for  X  with  respect  to  u,  one  can  easily  derive  that  the  turnover 

5 

angular  frequency  is  given  by 


u 


max 


R  C 
P  P 


V 


1  +  K 


(B7) 
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